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ABSTRACT 

Three series of rocket soundings were conducted from high latitude 
sites during winter. In the first series, four pitot pressure soundings 
were launched during a 13 hour period from Ft. Churchill, Canada 
(59°N) on January 31, and Feb. 1, 1967. The second series con- 
sisted of one pitot sounding and two acoustic grenade soundings 
carried out during a three hour period on January 13-14, 1970. 
Temperature and wind profiles and one density profile were 
observed independently to obtain the thermodynamic structure, the 
wind structure, and thus their interdependence in the mesosphere. The 
third series of soundings was conducted from Point Barrow, Alaska 
(71°N) on December 6, 1971. This series consisted of five sound- 
ings of which the first two and the last two were pitot-grenade pairs. 
Temperature profiles from all soundings in each series were aver- 
aged, and a smooth curve (or series of smooth curves) drawn through 
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the points. A hydrostatic atmosphere based on the average, meas- 
ured temperature profile was computed, and deviations from the 
mean atmosphere were analyzed in terms of gravity wave theory. 
The vortical wavelengths of the deviations were 10-20 km, and the 
wave amplitudes slowly increased with height. The experimental 
data were matched by calculated gravity waves having j period of 
15-20 minutes and a horizontal wavelength of 60-80 k7n. Our 
interpretation is generally consistent with the results of others 
who have studied gravity-acoustic waves in the atmosphere. The 
wind measurements are consistent with the thermodynamic meas- 
urements. The results also suggest that gravity waves travel 
from East to West with a horizontal phase velocity of approxi- 
mately 60 m sec“ l . 
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SYMBOLS 

to ~ circular frequency of the wave 
a = local sound speed 
k 0 - to/ a 

7 « ratio of specific heats 

g/a yjy- 1 = Brunt - Vaisala frequency 
W A~ 7g/2a~ a/2II 
g = acceleration due to gravity 
H = Pressure scale height 
A - to A / to 
B = co fl / to 
S = k x /lc 0 
G = cofo/to 
to h = a 2 /2k 

k = thermal diffusivity 
U = wind velocity of background field 
k = propagation vector 
k 2 = k 2 + k 2 

A it 

P = static pressure 
P° = stagnation pressure 
p - density 
u = velocity 
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INTRODUCTION 


Waves in a stratified fluid under the influence of gravity appear to have 
been initially discussed by Burnside (1889) and Love (1891). Both authors 
treated an incompressible fluid. Gortler (1943) used schlierin photography to 
show experimentally that disturbances in an incompressible stratified medium 
under the influence of gravity propagate along characteristic rays. Lamb (1909) 
treated a compressible, adiabatic, ideal gas whose density is stratified by 
gravity. Since the early contributions by Burnside, Love, and Lamb, there have 
been many theoretical papers on various types of gravity waves. Reasonably 
up-to-date treatments of the subject may be found in Eckart (1960) and Yih 
(1905). However, gravity wave theory is currently undergoing an active phase 
of development and contributions are being added to the literature at a rapid 
rate. A recent addition is the bool" by Tolstoy (1973). 

The present paper treats several series of rocket soundings of atmos- 
pheric temperature, density, and wind structure. The experimental results 
are interpreted in terms of gravity-acoustic wave theory. The experiments 
are of two types, pitot and grenade. For the pitot experiment, a pitot stagna- 
tion probe is mounted on the rocket vehicle. The position and velocity of the 
vehicle are measured by radar tracking. For an incompressible fluid, we 


have the simple formula 


where P° is stagnation pressure 
P is static pressure 
p Is density 
u is velocity 

Eq. (1) may also be written in the form; 


„ 2(P° - P L 
9 ~ u z 


(la) 


It is obvious from Eq. la that density is determined if stagnation pressure, 
static pressure, and velocity are measured. The addition of compressibility 
makes the analysis somewhat more involved without affecting the basic prin- 
ciple presented above. The pitot technique is described in detail by Horvath 
et. al (1962). 

For the grenade experiment, a series of grenades are exploded in the 
atmosphere. The position and time of the explosion is determined and the time, 
and the direction of arrival of the spherical sound pulse are measured by aground 
based microphone array. This information is used to deduce the mean temper- 
ature and horizontal wind speed in atmospheric slabs whose thickness varies 
between 2 km and 4 km depending on the particular grenade payload used. 
(Nordberg and Smith, 1964) 

Essentially, the pitot technique provides a high vertical resolution 
(0. 5 km) direct measurement of atmospheric density, while the grenade tech- 
nique gives a direct measure of atmospheric temperature and horizontal wind 
with a 2-4 km vertical resolution. For both techniques, the time to make the 
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atmospheric traverse is approximately one minute, resulting in a virtually 
instantaneous picture of the vortical structure. 

EXPERIMENTAL PROCEDURE 

The experimental data consist of three series of rocket soundings which 
were conducted from high latitude sites during winter. In the first series, four 
pitot soundings were launched during a thirteen hour period from Ft. Churchill, 
Canada (69°N) spanning January 31 and February l, 1967 (see Table 1), The 
data have been published by Smith, et al (1969). 

The second series consisted of three soundings carried out during n 
three hour period on January 13-14, 1970 from Churchill (Smith, et. al. , 1972) 
The first sounding was made with the acoustic grenude technique to measure 
the temperature and wind profiles in the 35-90 km region; the second sounding 
was a pitot probe launched 88 minutes later to measure the density profile; 
and the last sounding employed the acoustic grenade technique 104 minutes 
after the pilot sounding. Thus two temperature and wind profiles and one 
density profile were obtained independently to permit examination of the ther- 
modynamic structure, the wind structure, and the interdependence of each in 
the mesosphere. 

The third series of soundings was conducted from Point Barrow, Alaska 
(71°N) on December 6, 1971 (Smith et. al. , 1974). This series consisted of 
five soundings, of which the first two and last two were pitot-grenade pairs. 

The remaining sounding was performed with a pitot probe. The objective of 
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launching such pairs was to obtain high vortical resolution thermodynamic 
structure (pitot) and wind information (grcnndo) simultaneously. Tho sories 
commenced at 0300 GMT with a pitot probe, which was followed by a grenade 
at 0302 GMT. These profiles were flown over essentially the same trajectory 
only two minutes apart. The remaining soundings, a pitot at 0442 GMT, a pitot 
at 0752 GMT, and a grenade at 0802 GMT completed the series. 


OUTLINE OF GRAVITY WAVE THEORY 


Consider a series of layers of fluid, each lighter than the one below, 
as shown In Figure 1. A surface, wave in the lowest layer will transmit a 
disturbance to tho layer above it, which in turn will transmit the disturbance 
to tho layer above it, etc. Now let the density difference between layers and 
the layer thickness both approach small values e and 5, respectively. In the 
limit, as e and 8 approach zero, we have Internal gravity waves. The sound 
Ings described above have taken samples vertically through the layers. 

The simple theory involves linearizing the atmospheric equations of 
motion and specifying a sinusoidal solution, 1. e. , 

<p ~ <f> 0 exp P i( cot - k x x - k z z)l (2) 


Perturbations in pressure, density, temperature, and gravity wave 
generated wind are related to </> by what Hines (I960) has called polarization 
relations, and a dispersion relation which has the form: 






^ V • ••• 


S.'C. 


( 3 ) 
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where 


n x “ 


kx a 

CO 


n z = 


k a 
lv z 

co 


(4) 


(Pitteway and Hines, 1963) 

The following expressions aro obtained for phase velocity: 
co co <oX 2 _ X2 

Vpx = 5 V = isT c t 


(5) 


Also, 



provided that t > 


Since atmospheric gradients are usually gentle, the ray approximation 

for vertical wave propagation was used. Volland (1969a) has shown that the 

ray treatment is a sufficient approximation for gravity waves propagating 

obliquely upwards, provided the horizontal wave number is in the order of, or 

greater than co/a, where a is the sound speed, and co is the circular frequency 

co , 

of the ’’^ve, i. 0 . k x 3* — , or (co/a)/k x <1. Reference to eq. 5 shows that 
the last mentioned inequality may be written as V nv /a < 1. Thus the ratio 
of horizontal phase velocity to sound speed must be less than one. 

The above ratio, called VPXA, is plotted in Figure 2. We see that it 
is relatively insensitive to temperature and wave period, but varies linearly 
with vertical wave length. For vertical wave lengths less than 15 km, VPXA 
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is generally less than 0. 2. Thus, the condition for use of the ray approxima- 
tion is well satisfied. 

The next problem which comes to mind is how to take account of back- 
ground wind. The frequency of a moving fluid particle, or the 'intrinsic 
frequency' may be defined as 

to = w - k • U (6) 

where 

k is wave number 
U is wind velocity 

ce is wave frequency in a quiescent atmosphere 

The above value of intrinsic frequency may be substituted into equa- 
tions describing gravity wave propagation in a quiescent atmosphere, and a 
solution obtained. Jones (1969) discusses the conditions under whicn the above 
procedure is valid. Basically, the following assumptions must be satisfied: 

1. Vertical and horizontal derivatives of vertical winds must be less 
than (N 2 - ft 2 ) 1 / 2 , where N is the Brunt-Vaisala frequency. 

2. Horizontal derivatives of horizontal winds must be smaller than 

3. Time di ..Ivatives of horizontal winds must be much smaller than 
( to/ (3 ) x (N 2 - ft 2 ) 1 / 2 , where (3 = N 2 /g, 

4. Time derivatives of vertical winds are much smaller than g, the 
acceleration due to gravity. 
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Derivatives of winds are not generally known, but we can examine the 
range which is permitted by the above conditions. Condition "4" places a 
limit of 9.8 m sec" 2 to accelerations in vertical motion. Vertical winds 
generated by gravity waves are generally less than 10 m sec -1 . The fre- 
quency must be less than the Brunt-Vaisala frequency which is about 
0. 02 -1 . The resulting acceleration is 0. 2 m sec -2 , which is well below 
gravitational acceleration. Vertical wavelengths for major wind shifts are 
generally longer than 5 kilometers. Thus, we may approximate the maximum 
vertical derivative of vertical wind to be of the order of 10 m sec -1 per 5 km, 
or 2 m sec -1 km" 1 . Horizontal derivatives of vertical winds are less, since 
horizontal wavelengths are longer than vertical ones. For conditions con- 
sidered in this paper SI <0. 5N, consequently f2 z < 0. 25N 2 , and to a first 
approximation, we may simply require that vertical gradients of vertical wind 
be less than the Brunt-Vaisala frequency, i. e. 

dV /dz < N (7) 

Z 

Cancelling the space dimension, dV /dz becomes 0. 002 sec -1 . N is 

z 

approximately equal to 0. 02 sec -1 . Thus condition "1" is amply satisfied. 

Conditions "2" and "3" are the problematic ones. Taking £i 2 < N 2 , 
condition "3" simplifies to 
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Condition "2" may be written in the form 


dV 


H 


dx 


< ft 


(9) 


If time and space variations of horizontal winds are due to synoptic 
phenomena, then the following quantities are characteristic: 

V H ~ 60 m sec -1 ; t ~10,000 sec; dX ~500 km 

One can then estimate synoptic time and space derivatives to be 10 -3 m sec -1 , 
and 10“ 4 sec -1 . In the above case, eq. (8) requiros ( ft/N) to be greater than 
0. 001 and eq. (9) requires ft to be greater than 10 -4 sec -1 . 

ANALYSIS 03? DATA 

Table 1 shows the schedule of the four launches in the first series. 

The timos between launches are, respectively, 6h21m, 2h48m, and 3h32m. 
Temperature profiles from all soundings in each series were averaged and a 
smooth curve (or a series of smooth curves) drawn through the points. A 
hydro tatically determined atmosphere based on the average measured tem- 
perature profile was computed. The average temperature profile for the first 
set of data is given in Figure 3. Average background atmospheres were 
calculated in a similar manner for the second and third sets of experimental 


series. 


The pitot probe technique provides a measurement of the atmospheric 
density profile by relating the density to the measured ram pressure using 
appropriate aerodynamic theory, while the acoustic grenade technique provides 
a measurement of the atmospheric temperature profile by relating the tem- 
perature to the measured speed of sound using gas kinetic theory. The hydro- 
static equation and the equation of state are used to obtain the remaining two 
thermodynamic parameters in each technique. Such an approach assumes that 
vertical accelerations in the atmospheric motions are negligible relative to 
the acceleration due to gravity, and is valid for mean atmosphere calculations. 
Smith, et al (1968) compared the results of grenade soundings and pitot probe 
soundings at Wallops Island, Va. They found temperature agreement to be 
better than 3°K below GO km and better than 5°K between GO and 90 km altitude. 
Thus hydrostatic equilibrium was valid under the generally undisturbed con- 
ditions of Wallops Island, confirming the validity of the hydrostatic approx- 
imation for background atmosphere calculations. 

For the first test series, the ratio of the density measured by the pitot 
probe to the hydrostatically computed density (based on the average temperature 
profile) was plotted as a function of altitude for each probe. The density varia- 
tions are shown in Figure 4. The second and third soundings show a density 
minimum at around 110 km. This is most pronounced for the second sounding. 
The pronounced density drop, associated with a pronounced temperature rise, is 
quite probably due to an auroral event. The structure below 100 km contains 
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various features which could conceivably be traced from one sample to the next, 
but this typo of exercise will be left to the interested reader. Rather, the authors 
matched the experimental curves to results obtained from gravity wave theory. 
Our program was adapted from one used by Volland (1969 a, b). The periods which 
result in a reasonable fit are in a range of 15-20 minutes with a horizontal wave- 
length of 60-80 kin. The above values correspond approximately to those found by 
Gossard (19G2) who observed gravity waves in the troposphere. Figure 4 essen- 
tially gives only the vertical wavelength, \ z , thus only allowing one to obtain the 
product of cjX x . The reason foi Jiooslng the short period, short horizontal 
wavelength waves, is that these are more strongly damped. Reference to Fig- 
ure 4 shows that the growth of amplitude with altitude is extremely weak, while 
theory predicts that an undamped wave will grow exponentially with altitude. 

Figure 5 shows the spatial auto correlation function for atmospheric 
density between 28 and 80 km, and Figure 6 shows the corresponding power 
spectral density versus wavelength. Analysis techniques employed are dis- 
cussed in Bendat and Piersol (1966) and Blackman and Tukey (1958). (In this 
case "power" is defined in terms of wavelength and amplitude). Wavelengths 
between 2 km and 17 km are represented, with a primary peak around 14 km, 
and a secondary peak around 9 km. The auto correlation function for the second 
sounding shows a wavelength of 2x7km, and for the fourth, a wavelength of 2x9 
km, i. e. 14km and 18km, respectively, for the dominant sinusoidal component. 
Figure 7 shows representative gravity waves propagating in a sinusoidal wind 
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field, and in a quiescent atmosphere. The regular sinusoidal pattern present 
in the quiescent atmosphere density deviation with altitude is dramatically 
altered by the introduction of a regularly varying background wind. Thus in 
the real atmosphere, it is expected that the regular patterns superimposed on 
the thermodynamic structure by internal gravity waves will not be obvious. 

The interactions between propagating waves and varying background conditions 
are exceedingly complex and cannot be separated readily without a priori 
knowledge of the conditions. 

In the second test series, an attempt to identify propagating internal 
gravity waves In the mesosphere was made by sounding the region with three 
instrumented rockets on January 13-14, 1970, from Churchill, Canada 
(59°N). The three observations, which were conducted at approximately 
90 min. intervals, consisted of a pitot probe bracketed by two acoustic grenade 
soundings. Thus, two temperature and wind profiles and one density profile 
were obtained independently, permitting an examination of the thermodynamic 
structure and the wind structure. The first sounding was conducted at 2223 
GMT with the acoustic grenade technique to measure the temperature and 
wind profiles in the 35-90 km region; the second sounding was launched 88 
minutes later (2351 GMT) and employed a pitot probe to measure the density 
profile with 0. 5 km vertical resolution; and the last sounding again was made 
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with the acoustic grenade technique 104 minutes after the pitot sounding 
{0135 GMT). Thus, two temperature and wind profiles, and one density 
profile were obtained independently to permit an examination of the thermo- 
dynamic structure, the wind structure, and the interdependence of each in 
the mesosphere. For descriptions of the experiment, see Smith, et. al. 

(1972). 

An average atmosphere was calculated based on the three soundings, 
and the differences of the pitot data from the mean were computed (these are 
referred to as perturbations). The observed perturbations are shown as the 
solid curves in Figure 8 and are seen to give a pattern characteristic of in- 
ternal gravity waves in a stratified, compressible fluid. Recall that the ex- 
perimental density is measured directly, while the pressure is a dej ived 
quantity. In Figure 8, the period of the calculated gravity wave was 20 minutes, 
and the horizontal wavelength was 60 km. 

Figure 9 shows the comparison between the zonal wind profile meas- 
ured by the first grenade sounding, and the wind pattern associated with the 
theoretically calculated gravity wave. (A zonal drift wind varying linearly from 
12 m sec” 1 at 36 km to 44 m sec -1 at 96 km has been included). The wind pat- 
tern obtained by the second grenade sounding is generally similar. Figure 10 
shows the theoretically derived vertical velocity profile, which includes values 
of approximately 10 m sec -1 in the 80 km region. Justus and Edwards (1971) have 
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reported measured vortical velocities on the order of 20 m sec” 1 at altitudes 
between 88-118 km. Although these measurements wore made at a different 
latitude and time of year, thoy are at least comparable in magnitude. The 
Justus and Edwards data do confirm the existence of relatively largo vertical 
velocities at high altitudes. The vertical velocity profile given b.v Figure 10 
was included in the data reduction for the two grenade soundings and its effect 
on the mean profile was found to be small. 

The third set of measurements consisted of a series of pitot and grenade 
soundings from Point Barrow (71°N) launched on December G, 1971. The first 
two soundings, a pitot and a grenade, wore launched two minutes apart. The 
temperature perturbation obtali ed from the pitot profile is shown in Figure 11 
together with a theoretical gravity wave which was matched to it for a best fit 
(by eye). Figure 12 shows a comparison between the theoretical and observed 
wind patterns. The calculated wind pattern is that associated with the gravity 
wave whose temperature perturbation is given in Figure 11, A mean drift wind 
of 40 m sec' 1 has been included. The measured wind is obtained from the gre- 
nade sounding, so it represents a series of values in which the winds were aver- 
aged over layers (2-3 km thick) between grenade explosions. 

Figures 13 and 14 show the temperature perturbations derived from 
the observed density structure by the subsequent pitot soundings at 0442 GMT 
and 0752 GMT, respectively. The gravity wave structure is based upon the 
original match of Figure 11 and allowed to propagate with time to correspond to 
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the times of the observational data. Tho best match was found to have a 
horizontal wavelength of 70 1cm and a period of about 18 minutes. It is signif- 
icant >ihat only the first theoretical curve was fitted to the observed data, and 
that ;he remaining matches followed from the computer. Thu ft our data in- 
dicate that not only are the spatial density and wind structure measurements 
compatible with gravity wave theory, but also with the expected behavior with 
time. 

THE PARABOLIC TRAJECTORY EFFECT 

For simplicity of discussion, we have considered a vertical, instanta- 
neous slice through the atmosphere. In fact the actual sample is neither wholly 
vertical nor quite instantaneous, The latter approximation, however, is very 
good. The sampling time is one minute, representing only 5% of a wave period. 
The deviation of the actual trajectory from vortical is a more serious matter. 
Figure 15 shows a typical pitot probe trajectory. Figure 16 compares a 
trajectory going into the wave, a trajectory following the wave, and an instan- 
taneous vertical traverse. We see the expected Doppler phenomenon of short- 
ened wavelength going into the wave, and lengthened wavelength when going 
with the wave. 

The wave considered has a horizontal wavelength of 70 km and a period 
of 18. 5 minutes. While the measurements are being made, the horizontal 
distance traversed .by the rocket is 12 km, or l/6 of a wavelength. Thus, it is 
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the horizontal motion of the rocket rather than the finite sampling time which 
results in significant deviation from the original simplifying assumption. 

Examination of the experimental curves shows that the wavelength de- 
creases with altitude. Thus the rocket is going into tho wave train. Since the 
launch direction from Pt. Barrow is due east, we may conclude that tho 
wave train is coming from an easterly direction. The same situation pre- 
vails at Ft. Churchill, Our results suggest that tho waves in these instances 
travelled from East to West. The horizontal phase velocity is approximately 
60 m sec” 1 . 

CONCLUSIONS 

Three series of soundings were carried out at high latitude sites during 
winter. The first series gave four essentially instantaneous vertical density 
traverses during a 13 hour period at Ft. Churchill. The density variations from 
the mean stratospheric conditions were matched by gravity waves having ver- 
tical wavelengths of 14-17 km. A reasonable fit to the data was obtained with a 
plane gravity wave characterized by a horizontal wavelength of 70 km and a period 
of 18 minutes. The wave amplitude increased with altitude, but at a rate much 
slower than the exponential growth predicted for undamped waves. Thus one 
must conclude that the waves were heavily damped as they propagated upward 
through the atmosphere. 

Comparing our observations with the damping predictions of Pitteway 
and Hines (1963), it is found that for the constant viscosity case, (i. e. Increasing 
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kinematic viscosity) waves having the Input parameters considered, (i. d. 70 
km horizontal wavolength, and 10-20 km vertical wavelength) would tend to 
roach a peak amplitude at about 120 km altitude. The Pltteway-Ilines predic- 
tions are consistent with our observations, and their amplitude balance damp- 
ing case appears to be close to what we observed. The observations indicate 
a modest growth of amplitude with altitude, much closer to the no-growth 
approximation than to the exponential growth limit. 

Consider new the 70 km horizontal wavelength and the 18 minute period. 
Gossard (1962) obsorved gravity waves in the troposphere. He found the period 
to vary between 15 minutes and 120 minutes. The horizontal wavelength for 
waves of 15 minute period was 19 km, and for two hour period, it was 150 km. 
The difference between our observations and those of Gossard may be ex- 
plained by assuming that tropospheric waves with short vertical wavelengths 
are damped out. In fact Gossard (1962) found that the maximum leakage into 
the upper atmosphere occurs at periods of 10 to 20 minutes, starting with a 
white energy spectrum in the troposphere. Our value of 18 minutes is thus 
consistent with Gossard' s prediction.- Gossard found a horizontal phase speed 
of approximately 21m sec -1 . For an 18 minute period and 200°K, the follow- 
ing table relates our estimated horizontal phase velocity to vertical wavelength. 
200 degrees was chosen because it is approximately the temperature at the 


bottom of our test region. 
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20Q°K} 18 mlnuto period 


X y (km) 

X x (km) 

V px (m/sec) 

V pz (m/sec) 

5 

18.2 

16.9 

4. 6 

10 

36.4 

33.7 

9.3 

15 

54. 4 

60.4 

13.9 

20 

72.3 

67.0 

18.5 

25 

90.0 

83.3 

23.1 


From noctilucent cloud studies, Witt (1962) observed characteristic 
wavelengths of 30-40 km, and of 75 km, the latter being very nearly the 70 km 
deduced in the present study. Witt (1962) also estimated the wave velocity 
with respect to a frame of reference moving with the cloud system, obtaining 
a range of 70 to 135 m sec" 1 . The agreement with the values given in our tnblo 
above ifc reasonable. Noctilucent clouds are formed in the 80 to 85 Ian region 
of the atmosphere. 

Our results are also generally consistent with the discussion of internal 
atmospheric gravity waves at ionospheric heights presented by Hines (1960). 
Thus the 12 km vertical wavelength which Hines deduced from the meteor wind 
data of Greenhow and Neufeld (1955, 1959) is in reasonable agreement with our 
observations. Hines considered 125 m sec" 1 to be a typical horizontal phase 
speed and quoted the following experimental results. Munro (1958) found phase 
speeds to generally be in the range of 52-175 m sec" 1 , in agreement with both 
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our deductions and Wilt's (19 G2) noctilucent cloud observations, Iloisler (1958) 
found phase speeds to range between 97-207 m sec" 1 . The phase speed agree- 
ment is remarkable considering that ionosphoric heights range up to more than 
100 km above the top of our observation regime. Our period is more than ten 
times shorter than the 200 minutes deduced by Hines (1960) from correlation 
studies of meteor wind trails. However, a careful reading of Hines' paper 
suggests that the 200 minutes is an upper limit to the period* rather than a 
most characteristic value. 

The second series of soundingB is also compatible with a gravity wave 
interpretation. Vertical wavelengths were consistent with the values obtained 
in the first series, There was also reasonable agreement between measured 
and predicted horizontal velocities. 

The results of the third series confirmed the conclusions about vertical 
wavelength and amplitude of the first two series and the agreement of measured 
wind with gravity wave predictions. It was also possible to trace one gravity 
wave pattern through three consecutive samples, indicating that our period is 
correct. There is, of course, the possibility of aliasing. However, the slow 
wave growth with altitude points toward the short wavelength, short period 
solution, as do Witt’s observations at noctilucent cloud altitudes. 
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Table I 


Launch Sequence of First Series 
Fort Churchill, (59°N); 1967 


GMT 

Date 

Time After 
First Shot 

Type 

2317 

Jan 31 

OhOOm 

Pitot 

0538 

Feb 1 

6h21m 

Pitot 

0826 

Feb 1 

9h09m 

Pitot 

1158 

Feb 1 

12h41m 

Pitot 


Table II 

Launch Sequence of Second Series 
Fort Churchill (59°N); 1970 


GMT 

Date 

Time. After 
First Shot 

Type 

2223 

Jan 13 

OhOOm 

Grenade 

2351 

Jan 13 

lh28m 

Pitot 

0135 

Jan 14 

3hl2m 

Grenade 
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Table IH 


Launch Sequence of Third Series 
Point Barrow (71°N); 1971 


GMT 

Date 

Time After 
First Shot 

Type 

0300 

Dec 6 

OhOOm 

Pitot 

0302 

Dec 6 

0h02m 

Grenade 

0442 

Deo fi 

lh42m 

Pitot 

0752 

Dec 6 

4h52m 

Pitot 

0802 

Dec 6 

5h02m 

Grenade 




Illustration of Gravity Wave Concept and Nature of Experimental Measurement 




VPXA VPXA VPXA 


HORIZONTAL PHASE SPEED 





Figure 2. Gravity Wave Phase Speed as Function of Temperature, 
Vertical Wave Length, and Period 
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DENSITY DEVIATIONS, 1907 SERIES 
FORT CHURCHILL 



Figure 4 . Measured Density Structure for 1967 Ft. Churchill Series 
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CORRELATION 
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Figure 5. Spatial Auto-Correlation Functions for Density, 1967 Series 
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POWER DENSITY 
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Figure 6. Power Spectral Density Versus WavelenKth for 
Experimental Deviations, 19f>7 Series 
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Figure 7. Gravity Wave Comparable to Measured 
Vertical density Deviation for 1967 Series 
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ALTITUDE, KM 



Figure 8. The Observed Pressure and Density Perturbations (Solid Curves) 
Over Churchill on January 13-14, 1970 Compared With ifcft Perturbations 
Calculated From Theory (Broken Curves) as a Function t ,\ f Altitude 






ZONAL WIND. M SEC** 1 


Figure 9. The Observed Zonal Wind Component Compared 
With the Zonal Wind Profile Calculated From the 
Theory for the Described Experiment 
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TEMPERATURE PERTURBATION, °C 


Figure 11. 1971 Series Gravity Wave Matched 
to First Pitot Temperature Deviation 






TITUDE, KM 



Figure 12. Comparison Between Theoretical and Observed 
Wind Patterns for 1971 Series 
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ALTITUDE, KM 



TEMPERATURE PERTURBATION, °C 


Figure 13. Comparison of Second Pitot Temperature Deviation 
with Gravity Wave Prediction for 1971 Series 


a 
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TEMPERATURE PERTURBATION, °C 

Figure 14. Comparison of Third Pitot Temperature Deviation 
with Gravity Wave Prediction for 19 G7 Series 
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ALTITUDE ABOVE 33 km 



PARABOLIC TRAJECTORY EFFECT ON GRAVITY WAVE 



Figure 16. Gravity Wave Doppler Effect 



